function rate=massFunction(mass)

%hurley 2008
m0=0.1;
if(mass>1)
    rate=mass^(-2.3);
else
    rate=NaN;
end

end